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ABSTRACT 


Heat transfer, thermal stresses, and thermal buckling analyses were performed on the unconventional 
wing structures of a Hyper-X hypersonic flight research vehicle (designated as X-43) subjected to 
nominal Mach 7 aerodynamic heating. A wing midspan cross section was selected for the heat transfer 
and thermal stress analyses. Thermal buckling analysis was performed on three regions of the wing skin 
(lower or upper); 1) a fore wing panel, 2) an aft wing panel, and 3) a unit panel at the middle of the aft 
wing panel. A fourth thermal buckling analysis was performed on a midspan wing segment. The unit 
panel region is identified as the potential thermal buckling initiation zone. Therefore, thermal buckling 
analysis of the Hyper-X wing panels could be reduced to the thermal buckling analysis of that unit panel. 
“Buckling temperature magnification factors” were established. Structural temperature -time histories are 
presented. The results show that the concerns of shear failure at wing and spar welded sites, and of 
thermal buckling of Hyper-X wing panels, may not arise under Mach 7 conditions. 

NOMENCLATURE 

a length of rectangular plate, in. 

b width of rectangular plate, in. 

CL clamped boundary condition 

c p specific heat, Btu/lb- 0 F 

E Young's modulus, lb/in 2 

E23 two-nodes bar element 

E33 triangular combined membrane and bending element 

E43 quadrilateral combined membrane and bending element 
FR panel edge supports allowing free in-plane motion (free) 

FX panel edge supports with no in-plane motion (fixed) 

JLOC joint location (node) 

k thermal conductivity, Btu/hr-ft- 0 F 

LTA Lockheed Thermal Analyzer 

SPAR Structured Performance and Resizing 

SS simply supported boundary condition 

T temperature, °F 

T mat average temperature for material properties, °F 

T n temperature at node n 

T a peak temperature of dome-shaped temperature distribution, °F 

T s panel boundary heat sink temperature, °F 

T u uniform temperature distribution, °F 
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Dryden computer simulation program 
thermal protection system 

time counted from the air-launching of Pegasus booster rocket from B-52, sec. 

computer program to compute view factors, see equation (1) 

critical buckling value at buckling 

coefficient of thermal expansion, in/in- °F 

emissivity 

= ( T o ) _/( T u ) ^Buckling temperature magnification factor of the first kind = buckling 
temperature (T 0 ) cr of dome- shaped temperature loading (fixed or free support condition) 
divided by the buckling temperature (T u ) cr of uniform temperature loading (fixed 
support condition). 

eigenvalue associated with uniform temperature loading 
eigenvalue associated with dome-shaped temperature loading 
Poisson's ratio 

= 0 ) cr^ (T 0 ) cr I T = 0 Buckling temperature magnification factor of the second 

kind = buckling temperature (T 0 ) cr of dome- shaped temperature loading case with any heat 
sink temperature T s * 0 (free support condition) divided by the buckling temperature (T 0 ) cr 
of dome- shaped temperature loading case with zero heat sink temperature T s = 0 (free 
support condition) 
density, lb/in 3 

uniaxial tensile failure stress, lb/m" 
ultimate tensile strength, lb/in 2 
chord-wise stress, lb/in 2 
yield strength, lb/in 2 
sheer failure stress, lb/in" 
maximum shear stress, lb/in" 


INTRODUCTION 

Hypersonic flight vehicles are subjected to severe aerodynamic heating during flights. To maintain 
structural integrity at these high temperatures, the vehicle structural design concepts for high Mach 
number vehicles are different from those of low Mach number aircraft. The vehicle structure may be 
called “hot” structures or “warm” structures, depending on the operating temperature range. The hot 
structures are fabricated with high-temperature materials and are capable of operating at elevated 
temperatures exceeding 1000 °F. Typical hot structural components are carbon/carbon composite 
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structures, hat- stiffened panels fabricated with either monolithic titanium alloys or metal-matrix 
composite materials, and honeycomb sandwich panels fabricated with high-temperature materials such as 
Titanium, or nickel-based Inconel® (Inco Alloys International, Inc., Huntington, West Virginia) alloys. 
The warm structures are fabricated with light-weight materials such as aluminum and must be insulated 
so that the sub-structural temperatures will not exceed the operating temperature limit of 350 °F. The 
space shuttle orbiter is a good example of warm structure. The entire vehicle is protected with a thermal 
protection system (TPS) to shield the aluminum substructure from overheating beyond the warm 
temperature limit (ref. 1). 

An example of a recent hot structure is the new hypersonic flight research vehicle called Hyper-X 
(designated as the X-43 vehicle), which has unconventional wing structures with irregular-shaped wing 
panels (described in the following section). 

Temperature distributions in the hot structures are, in general, not uniform because of cooler 
substructures such as spars and ribs. The nonuniformity of the structural temperature distribution induces 
thermal stresses. Excess thermal stresses could cause vehicle wing panels to buckle and creep, and could 
also cause the skin and spar weld sites to lose bonding integrity. A uniform temperature field induces no 
thermal stresses in structures fabricated with materials that have the same coefficient of thermal 
expansion. Thermal buckling problems are of great concern for both hot and warm structures, and have 
been extensively studied by Ko (refs. 2-11). Those past thermal buckling studies were conducted on 
rectangular- shaped honeycomb sandwich and hat-stiffened panels, and not on irregular- shaped panels 
such as the Hyper-X wing panels. 

To fully understand the thermo structural performance of the unconventional Hyper-X wing structures 
under nominal Mach 7 aerodynamic heating , heat transfer analysis, thermal stress analysis, and thermal 
buckling analysis must be performed. This requires complex three-dimensional modeling of the Hyper-X 
wing structures. However, without going through complex modeling of the entire wing structure, one can 
use a simple method to quickly approximate the panel buckling temperatures of wing panels. This report 
presents the results of heat transfer analysis and thermal buckling analysis of wing structures using such 
an approximation method, and discusses how the Hyper-X wing structures would perform under nominal 
Mach 7 heating. 

Use of trade names or names of manufacturers in this document does not constitute an official 
endorsement of such products or manufacturers, either expressed or implied, by the National Aeronautics 
and Space Administration. 

HYPER-X FLIGHT RESEARCH VEHICLE 


Hyper-X (designated as X-43) is a new hypersonic flight research vehicle (12 ft long, 5 ft span, 
3,000 lb weight), designed to be flown at a range of Mach 7~10 (fig. 1). The underside of the vehicle 
consists of an inlet ramp, a scramjet (supersonic-combustion ramjet) engine module, and an expansion 
ramp. The scramjet engine is an air-breathing ramjet engine in which the airflow through the whole 
engine remains supersonic. The Hyper-X scramjet engine uses hydrogen as fuel and inlet air as oxidizer 
(no oxidizer is carried on board). Tire proposed flight trajectory of Hyper-X is shown in figure 2. Tire 
Hyper-X rides on a winged Pegasus booster rocket, which is carried under the wing of a B-52 aircraft up 
% 

It should be noted that the maximum Mach number reached during the nominal Mach 7 mission was 7.5. 
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to a launch altitude of 17,000 ft. for the Mach 7 mission or 43,000 ft for the Mach 10 mission. After air 
launching from the B-52, the Pegasus booster rocket will accelerate and ascend to an altitude of 
approximately 100,000 ft, reaching the test velocity (of Mach 7 or Mach 10). After separation from the 
booster rocket, the cowl door of the Hyper-X scramjet engine opens to test the performance of the 
scramjet engine. Once the cowl door is open, fuel is injected, ignited and burned for about 8 seconds. The 
entire event from the opening to the closing of the cowl door lasts for 34 seconds. 

The Hyper-X wing structure is fabricated with high-temperature Haynes 230 alloy (a 
nickel-chromium-tungsten-molybdenum alloy) which has relatively low thermal expansion 
characteristics. The design concept of its wing structures (fig. 3) is entirely different from that of the 
conventional wing structures. The conventional spar and rib system is replaced with multiple radial 
stiffeners (spars, 0.25 in. wide) fanning out from the pivoting wing roots. To house the instrumentation 
inside the wing structure, upper and lower wing skins (0.090 inches thick) are divided into two separate 
wing panels (a fore wing panel and an aft wing panel). The wing panels are then butt- welded at their 
edges to the main whig frame, and line-welded to the radial reinforcing stiffeners without using 
conventional fastening screws or rivets. Because the edges of the heated wing panels are constrained, 
potential thermal buckling of the wing panels and possible shearing off of the line-welded sites are of 
great concern. 


DESCRIPTION OF PROBLEM 


Heat transfer and thermal stress analyses are to be performed on the midspan segment of the Hyper-X 
wing (fig. 3). A nominal Mach 7 aerodynamic heating is to be used as heat input to calculate time histories 
of structural temperatures. The worst-case non-uniform temperature distribution obtained from the heat 
transfer analysis is to be used for temperature loads to calculate thermal stresses in the wing segment. 

Thermal buckling analysis is to be perfonned on the following three regions of the wing skin panels 
(lower or upper): 1) the fore wing panel, 2) the aft wing panel, and 3) a unit panel at the middle of the aft 
wing panel (fig. 3). hr addition, thermal buckling is also to be conducted for the wing segment (fig 3). 
These analyses are being done to locate the thermal buckling initiation zone. 

The unit panel is used to establish buckling temperature magnification factors (1) to scale up the 
buckling solution of uniform temperature loading for generating a buckling solution for dome- shaped 
temperature loading; or (2) to scale up the buckling solution of dome-shaped temperature loading with 
unheated boundary heat sinks for generating a buckling solution for the same case when the heat sinks are 
heated up. 


HEAT TRANSFER ANALYSIS 


The finite-difference heat transfer analyses of Hyper-X wing structures are perfonned using the 
Lockheed Thenual Analyzer (LTA). The LTA is a generalized FORTRAN program that computes 
transient temperature distributions in complex structures based on an electrical analogy of capacitors and 
resistors. 
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Finite-Difference Thermal Modeling 


The Hyper-X wing segment selected for the heat transfer analysis is located at the wing midspan, 
approximately 4.7 in. from the wing root edge (fig. 3). The chord-wise region lies between the 
stream- wise distances 11.0 in. and 18.0 in. measured from the carbon/carbon leading edge, and spans 
over three neighboring radial stiffeners (or spars) (fig. 3). 

Two finite-difference thermal models (fig. 4) were generated for two cases of welded site skin and 
spar contacts; 1) full contact, and 2) partial contact (inset of fig. 4). The two thermal models are identical 
(fig. 4) except for the contact resistance between the skin panels and radial stiffeners. Each model 
consists of 64 capacitors modeling Haynes® 230® (Haynes International, Kokomo, Indiana) alloy. Outer 
structural mold lines use a 1.5 degree half- angle for both upper and lower skins. For the full contact, the 
wing skin panels are perfectly bonded to the full width (0.25 in.) of the stiffeners. For the partial contact, 
the welded site is 0.040 in. wide and 0.01 in. thick. The resistance between capacitors { 1 and 66}, {13 
and 70}, {25 and 75}, {69 and 41}, {74 and 53}, and {79 and 65} (fig. 4) were calculated for 1) full 
contact conductance, or 2) partial contact conductance representing a 0.040-inch- width contact area. The 
partial contact case is a more realistic representation of the actual skin and spar welded sites. The thermal 
models have a surface emissivity of e= 0.85. The external radiation view factors were calculated by 
hand, but the internal radiation view factors were calculated from the VIEW program (described below). 
The input thermal properties of Haynes 230 alloy are listed in Table A-l of the Appendix. 

Internal Radiation View Factors 

The internal radiation view factors were calculated from the VIEW program, which is incorporated 
into the Structural Performance And Resizing (SPAR) finite-element computer program (ref. 12). Figure 
5 shows the two SPAR finite-element radiation models generated for the internal view factor 
computations. Each model has 63 internal radiation elements. Each of the radiation elements matches the 
interior surface of each corresponding element of the finite-difference thermal model (fig. 4). 

In the computations of internal radiation view factors using the VIEW program, two cases were used: 
the combined case and the separate case. For the combined case, the entire wing segment was considered 
as a single enclosure with obstruction at the middle spar (fig. 5-a). For the separate case, each bay was 
considered independently as an enclosure without obstructions (fig. 5-b). 

For case one, the radiation element numbering was stalled from bay 1 and continued on bay 2, (not 
independently numbered for each bay). The command OBSTRUCT was used to define the obstructing 
surfaces, which are the two sides of the middle spar (i. e., elements 12, 13, 14, 15, 16, 59, 60, 61, 62, 63) 
(fig. 5-a). The ENCLOSURE command was used to identify the system of radiation surfaces as an 
enclosure, and to correct the calculated view factors so that the sum of the view factors from each 
radiation element is equal to unity. Namely, if F fJ (/, j = 1 , 2, 3, . . ..) is the radiation view factor defined as 
the fraction of radiant heat, leaving radiation element i and incident on radiation element j, then the 
ENCLOSURE command will enforce the following condition for the final values of F ;/ . 

ZO = 1 (i) 

j 

where the summation is taken over j for each given i. 
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For case two, the internal radiation view factors were computed independently for each bay using the 
VIEW command. For this case, the radiation elements are numbered independently (stalling from No. 1) 
for each bay (fig. 5-b). Because there is no obstruction within each bay, the ENCLOSURE and 
UNOBSTRUCT commands were used. Tire command OBSTRUCT was not needed. The above two 
approaches (fig. 5-a, 5-b) gave identical internal view factor solutions, thus verifying the accuracy of 
view factor calculations. 


Heating Profile 

Figure 6 shows the preliminary flight trajectory for tire Mach 7 mission of the Pegasus booster rocket 
carrying the Hyper-X vehicle. The maximum Mach number of the flight trajectory turned out to be 7.5 
(fig.6-c). This was the flight trajectory available at the time of this analysis and was, therefore, used to 
calculate the aerodynamic heating rates. Tire Dryden in-house computer code called TPATH was used for 
these calculations. Various parameters for inputs to the TPATH code are: time histories of altitude, angle 
of attack, and Mach number (fig. 6) as well as the outer mold-line geometry of the wing cross section. 
The program calculates transient heating rates and surface temperatures, and also computes heat transfer 
coefficients, boundary layer recovery temperatures and other parameters required to calculate the 
aerodynamic heating rates. The program permits the use of different theories for calculating the heat 
transfer coefficients. These theories can be properly applied at each location of interest for laminar or 
turbulent flow condition, or for flows with transition. The transition condition can be input as a function 
of Reynolds number and local Mach number, or prescribed at a specific time. In the present analysis, the 
transition criteria used resulted in turbulent flow calculations for the entire flight profile. Local flow 
conditions were calculated for an attached flow using the oblique shock theory (ref. 13). Heat transfer 
coefficients were calculated using Eckert’s reference enthalphy method (refs. 14, 15). Calculations were 
made for the upper and lower surface of the horizontal wing. Real gas properties of air were used in all 
calculations (ref. 16). Heat transfer coefficients and recovery temperatures calculated from the TPATH 
code were then used as inputs to the LTA. This thermal analyzer program then calculated the heating 
rates and corresponding structural temperatures for the two-dimensional thermal models. 

THERMAL STRESS ANALYSIS 


The SPAR finite-element computer program (ref. 12) is used in the linear thermal stress analysis of 
the wing segment. The following sections discuss finite element modeling, boundary constraints, and 
temperature load. 


Finite Element Modeling 

The wing segment selected for the thermal stress and thermal buckling analysis is located at the whig 
midspan cross section where the heat transfer analysis is performed (fig. 3). The wing segment has a unit 
span-wise width, and spans over three neighboring radial stiffeners (spars) in stream- wise direction. The 
finite-element structural model generated for the wing segment is shown in figure 7. The wing panels and 
the spars are modeled with a single row of quadrilateral combined membrane and bending elements (E43 
elements). The nodal coordinates of the finite-element structural model are made coincidental with those 
of the finite-difference thermal model (fig. 4). Thus, the nodal temperature output from the thermal model 
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can be used directly as temperature input to the structural model. Only one welded site at the top of the 
middle spar (radial stiffener) is modeled in detail for shear stress calculations (inset of figure 7). The wing 
skin and spar welded site is modeled with E43 elements whose lower boundaries are connected to the 
middle spar (E43) element through triangular membrane and bending elements (E33 elements). To 
simulate the effect of stream-wise thermal expansion restraint resulting from the cooler wing frame, the 
nodes of the three radial stiffeners are interconnected with two nodes (E23 bar elements). The size of the 
finite-element structural model for the wing segment is listed in Table 1. 


Table 1. Number of joint locations and elements of 
wing- segment finite-element model. 

Item 

Number 

JLOC 

166 

E23 elements 

20 

E33 elements 

8 

E43 elements 

95 


The temperature-dependent material properties of Haynes 230 used for input to the structural model 
are listed in the Appendix. 


Boundary Constraints 

The wing-segment model is allowed to expand freely in all directions. The span-wise edges of the 
wing panels are constrained to have zero slopes in the span-wise direction to simulate the effect of 
continuity to neighboring regions. 


Temperature Load 

The structural temperature distribution at time t = 89 sec (from launch of the Pegasus booster) 
obtained from the heat transfer analysis is used as thermal load input for the thermal stress analysis. Time 
r = 89 sec is the instant when the difference between the upper and lower skin peak temperatures reaches 
a maximum. As will be seen later, this structural temperature distribution over each bay is arch- shaped. 
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THERMAL BUCKLING ANALYSIS 


The SPAR program is also used for the linear thermal buckling analysis of the wing panels. The 
following sections discuss material property iterations, finite element modeling, boundary conditions, 
temperature load, and buckling temperature. 

Material Property Iterations 

In thermal buckling analysis, when high accuracy buckling solutions are sought, the material 
property iteration process is required. This iteration process is described in great detail in ref. 17. 
Initially, the buckling temperature is still unknown and, therefore, the material properties at a certain 
assumed temperature (usually room temperature) are first used as input. Once the first-iteration buckling 
temperature is calculated, the material properties corresponding to the first calculated buckling 
temperature are input to the second iteration of eigenvalue calculations. This material iteration process 
(updating material properties) must continue until the temperature used for the input material properties 
converges with the final buckling temperature. In the present buckling study, the material iteration 
process (averaging 5 ~ 7 iterations) for each case was terminated when the assumed material temperature 
approached the calculated buckling temperature with less than 0.5 percent error. 

Linite Element Modeling 

The wing-segment model (fig. 7) is also used for thermal buckling analysis. Figure 8 shows the three 
finite-element models (called the fore panel model, the aft panel model, and the unit panel model ) 
generated respectively for the fore wing panel, aft wing panel, and the unit panel located at the middle of 
the aft wing panel (fig. 3). All the wing panels are modeled with E43 elements. The boundary heat sinks 
and radial stiffeners are modeled with E23 bar elements. The sizes of these three wing finite-element 
models are listed below. 


Table 2. Number of joint locations and elements of wing-panel 
finite-element models. 



Fore panel model 

Aft panel model 

Unit panel model 

JLOC 

3116 

2419 

984 

E23 elements 

350 

276 

126 

E43 elements 

3000 

2320 

920 
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Boundary Conditions 


Different support conditions listed below are considered for comparative studies of how the buckling 
temperatures change with the support conditions. 

Table 3. Definition of support conditions. 



Panel Boundaries 

Panel and Stiffener Welded Sites 

SS-SS 

simply supported 

simply supported 

SS-CL 

simply supported 

clamped 

CL-SS 

clamped 

simply supported 

CL-CL 

clamped 

clamped 


Based on the fact that the Hyper-X wing panels (fig. 4-b) are line-welded, the SS-SS condition listed 
in Table 3 could be the closest to the actual support condition of the Hyper-X whig panels. 

Temperature Load 

The thermal load inputs used are uniform temperature loading and dome- shaped (or arch -shaped ) 
temperature loading. 

Uniform Temperature Loading 

Because the actual temperature distributions over the entire wing surfaces are not available, uniform 
temperature loading will be analyzed first. Figure 9(a) shows an example of a typical uniform 
temperature loading case for the unit panel. The uniform nodal temperature input is chosen to be 
T = T = 1 °F . The eigenvalue X c (scaling factor) calculated from the SPAR program will then give the 
buckling temperature (T u ) cr Namely, 


(T , ) = A x 1°F (2) 

This buckling solution is for the uniform temperature loading over the wing panels with fixed 
supports, and is not the actual case. This buckling solution is the simplest fundamental case, and can 
show the location of the weakest bay, where thermal buckling is likely to take place. When the panel 
boundaries can have free in-plane motion (free supports), the uniform temperature loading case obviously 
can not induce thermal buckling. 

Dome-Shaped Temperature Loading 

Hot structural panels are usually fastened to the cooler substructures (of relatively large mass) that 
function as heat sinks. Thus, even under a uniform surface heat flux, the temperature distribution over the 
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hot panels will not be uniform, but dome-shaped (ref. 18). This is the typical behavior of hot structural 
panels. 

For the whig segment, the calculated dome-shaped temperature distribution obtained from the 
wing-segment heat transfer analysis is used as the temperature load input. 

For the two-dimensional wing panels, temperature distribution over the wing surface must be 
calculated from heat transfer analysis by the modeling of one complete wing structure, hi order to avoid 
that complexity, the following quick approximation method will be used. 

Because the cooler spars function as heat sinks, the calculated temperature distribution over each bay 
is expected to be dome-shaped based on past hot structural experiments (ref. 18). By observation (fig. 3), 
the unit panel that has the largest area compared with other sub-panels must be the weakest sub-panel hi 
resisting thermal buckling. Therefore, buckling analysis of the entire wing panels under actual 
temperature loading can then be reduced to the analysis of the unit panel that is under dome- shaped 
temperature loading to obtain approximate buckling solutions. 

As will be seen shortly, the stream-wise distribution of panel temperature over each bay exhibits an 
arch- shaped curve (a cross sectional shape of a dome surface) because of the existence of heat sinks 
(cooler radial stiffeners). Such an arch-shaped curve is very similar to part of the sine curve. Therefore, 
for the two-dimensional case, dome-shaped temperature distribution could be a good approximation of 
the actual surface temperature distribution over any unit panel that is supported by boundary heat sinks. 

For a rectangular plate with length a and width b, the dome-shaped temperature distribution (with 
peak temperature T 0 and boundary heat sink temperature T s ) may be ideally described by the following 
equation. 


T(x,y) = r, + <r o -r s )siAinS 


(3) 


which gives a sine x sine surface lifted by an amount T s . 

Equation (3) may be applied directly to the irregular-shaped unit panel by simply distorting the 
boundaries and mesh of the finite-element model generated for the rectangular plate to form the irregular 
boundaries and mesh of the unit panel model. Thus, the nodal temperatures T(x,y) calculated from 
equation (3), will automatically provide the nodal temperatures for the unit panel model. Figure 9(b) 
shows the dome-shaped temperature loading over the unit panel calculated from equation (3). For the 
Hyper-X case, the ratio of averaged heat sink temperature T s to the wing panel maximum temperature T a 
is roughly T S /T o = 0.54 at tune t = 89 sec. 

For the dome-shaped temperature load input to the SPAR program for the eigenvalue calculations, the 
peak temperature T o = 1 °F was used. The heat sink temperature T s was allowed to vary over the range 
T/T 0 = 0~1 (i.e„ T/T o = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.54, 0.6, 0.7, 0.8, 0.9, 1.0; with 
T S /T o = 0.54 being the Hyper-X case). The buckling temperature (T 0 ) cr for the dome-shaped 
temperature loading case is then obtained from 

(T ) = l x 1°F (4) 

v o'er o v 7 
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where X () is the eigenvalue calculated from the SPAR program for the dome- shaped temperature loading 
case. The reason for using the whole range of T /T o = 0~ 1 is to study the effect of the heat sink 
temperature T s on the panel buckling temperature (T 0 ) cr Keep in mind that the panel buckling 
temperature (T 0 ) cr will increase with the increasing heat sink temperature T s as a result of relaxation of 
the thermal expansion constraint exerted on the unit panel by the boundary heat sink. 


Buckling Temperature Magnification Factors 


As will be seen shortly, the buckling temperature (T 0 ) cr for the dome-shaped temperature loading 
case is much higher than the buckling temperature (T u ) cr of the uniform temperature loading case (with 
fixed supports). Therefore, a buckling temperature magnification factor of the first kind rj defined below 


V 0 ) 


(5) 


will be used to indicate how many times the buckling temperature (T 0 ) a . of the dome- shaped temperature 
loading case (free or fixed supports) is magnified from the buckling temperature (T c ) cr of the uniform 
temperature loading case with fixed supports (fundamental case). 

In reality, the wing panel attempts to expand under heating, but its expansion is resisted by the cooler 
boundary substructures (heat sinks) which expand less. This boundary constraint is the cause of thermal 
buckling of the wing panels. Such substructure constraints will gradually relax as the substructures are 
heated up, resulting in a higher panel buckling temperature. In order to find out how the heat sink 
temperature T s affects the panel buckling temperature (T 0 ) cr for the dome- shaped temperature loading 
case with free supports, another buckling temperature magnification factor of the second kind q defined 
below 




(T ) 

v o'er 


( T o\r I T s = 0 


(6) 


will be used to indicate how many times the buckling temperature ( T 0 ) cr for the dome-shaped temperature 
loading case with free support (SS or CL) is magnified when the heat sink temperature T s increases from 
7' v = 0 (no heat sink thermal expansion) to a certain non-zero value T^O (with heat sink thermal 
expansion). 


RESULTS 


The results of the heat transfer analysis, thermal stress analysis, and thermal buckling analysis of the 
Hyper-X wing structure subjected to Mach 7.5 heating (the nominal Mach 7 mission) are presented in the 
following sections. 
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Heat Transfer 


Figure 10 shows the time histories of nodal temperatures { 7j 3 (t), T 10 (t)} at nodes {13, 70} of the 
middle spar upper welded site (fig. 4) for the full and partial contact cases. The nodal temperature 
difference (T 13 - T 10 ) is slightly higher for the partial contact case, and it reached a maximum value of 
(T i3 — T ?q) = 150°F at approximately t = 83 sec. into the flight profile. 

Figure 11 shows the time histories of nodal temperatures {T 53 (t), T 14 (t)} at nodes {53, 74} of the 
middle spar lower welded site (fig. 4) for the two cases of contact conditions. Although the nodal 
temperature difference (T 53 - T 14 ) at the lower weld site is less severe than (T 13 - T 10 ) of the upper weld 
site, it is still substantial. For the partial contact case, the lower welded site nodal temperature difference 
(T 53 - T 14 ) reached a maximum value of (T - 3 - T 74 ) = 130 °F at approximately t = 123 sec into the 
flight profile. 

Shown in figure 12 are temperature time histories at locations down the middle spar and along the 
upper skin panel. The temperature gradient along the skin panel (between capacitor 13 and 14) appeal’s to 
be greater than the temperature gradient from capacitor 13 to 70, but the distance between the two 
capacitors is longer. The temperature gradient per inch between capacitor 13 and 70 is actually larger. 
Temperature time histories on the upper skin panel can be seen in figure 13. The temperature gradient 
between capacitor 13 and 15 is higher than the gradients between capacitors 13 and 17 or 13 and 19. 
Distance between capacitors 13 and 15, 15 and 17, and 17 and 19 is 0.5 inches apart. Figure 14 shows 
nodal temperature time histories for node 7 of bay 1 upper skin and node 46 of bay 1 lower skin. Note 
that AT = Tj - T 46 reached a maximum value, AT = 546 °F at t = 89 sec. As can be seen in the next 
figure, node 7 and node 46 are respectively the peak temperatures of the upper and lower skins. 

Figure 15 shows the chord- wise distributions of the upper and lower skin temperatures at t = 89 sec 
into the flight profile. These temperature distributions were calculated from the thermal model with 
partial contact at the spar and skin weld sites. As can be seen from figure 15, the temperature distributions 
are steepest near the spars (heat sinks). Note that the panel temperature distribution over each bay is 
arch-shaped. This is the typical behavior of hot structural panels supported by the boundary heat sinks. 
The structural temperatures at typical points of the wing segment at t = 89 sec are listed in Table 4. The 
upper and lower wing skin temperatures listed are the peak temperatures of each bay. 


Table 4. Structural temperatures at typical points of wing segment; r = 89 sec, °F 



Fore 

Bay 1 skin 

Middle 

Bay 2 skin 

Aft 


weld site 

(peak) 

weld site 

(peak) 

weld site 

Upper 

675 

1,153 

748 

1,127 

638 

Lower 

479 

607 

476 

592 

418 


Notice that the designed flight trajectory caused the Hyper-X upper skin temperatures to be 
significantly higher than the lower skin temperatures. This implies that buckling is more critical for the 
upper wing panels than the lower wing panels. The panel temperature distributions of figure 15 and the 
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associated spar temperatures (not shown) were used as temperature inputs to the structural model for 
thermal stress and thermal buckling analyses of the wing segment. 

Thermal Stress Analysis 

The calculated chord-wise thermal stresses cr Y induced in the wing-segment skins under free 
expansion (to simulate actual situation) are listed in Table 5. 


Table 5. Chord-wise stresses o x in the wing-segment skins; free expansion, lb/im. 



Front bay 

Rear bay 

Upper skin 

-64,900 

-61,364 

Fower skin 

-37,419 

-31,938 


Because the chord- wise thermal expansion of the wing panels are restrained by the cooler wing frame 
(heat sink), both upper and lower wing panels of each bay are under compression. Even though the 
temperature distribution over the wing panel of each bay is non-uniform and arch-shaped (fig. 15), the 
thermal stress induced in the wing panel of each bay is constant. This is the typical behavior of hot 
structural panels. 

Figure 16 shows the deformed shape of the welded site under compression from wing panels of 
opposite bays, with element shear stresses indicated. The maximum shear stress T max is located at the 
aft upper corner element of the welded site (fig. 16) with the value T max = 38,564 lb/in 2 . The 
temperature at the middle spar upper welded site was calculated as T = 748 °F at t = 89 sec. The 
shear failure stress zy at this temperature level is approximately 1 y ~ 54,000 lb/in 2 for Haynes 230 
alloy. This value is one-half of the uniaxial tensile failure stress Of = 108,000 lb/in2 at T = 748 °F 
based on Mohr's circle. Based on this current simple stress analysis, the maximum shear stress T max at 
the weld site will reach as high as 71 percent of the shear failure stress zy during Mach 7.5 heating. 

Thermal Buckling Analysis 

The results of thermal buckling analysis of the Hyper-X wing panels are presented in the following 
sections. On average, about 5~7 material iterations were required to obtain acceptable buckling solutions 
presented in these following sections. 

Uniform Temperature Loading 

Figures 17, 18, and 19, respectively, show the buckled shapes of the fore panel, aft panel and unit 
panel under uniform temperature loading and under different fixed support conditions. Note that, except 
for the SS-CF condition, the aft panel buckled at the unit panel region, and the buckled mode in this 
region (fig. 18) is similar to the buckled mode of the isolated unit panel (fig. 19). The buckling 
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temperatures calculated for those wing panel models subjected to uniform temperature loading are shown 
in Table 6. 


Table 6. Buckling temperature of different wing-panel models 
under uniform temperature loading; fixed supports, (T u ) cn °F. 


Support 

condition 

Fore panel 

Aft panel 

Unit panel 

SS-SS 

167 

121 

97 

SS-CL 

176 

188 

239 

CL-SS 

355 

175 

114 

CL-CL 

410 

294 

295 


The buckling temperatures listed in Table 6 are for the uniform temperature loading over the whig 
panel and not for the actual temperature loading. Because of the fixed boundary supports and fixed panel 
and stiffener joint site, the buckling temperatures are amazingly low. Note that, except at the SS-CL 
condition, the aft panel has lower buckling strength than the fore panel, and the buckling critical region is 
located at the unit panel region of the aft panel (fig. 18). 

In reality, the panel support sites (heat sinks) can expand when heated, resulting in the reduction of 
thermal stress levels in the wing panels as a result of reduced heat sink restraint. Therefore, the actual 
buckling temperatures could be much higher than the values presented in Table 6. 

Dome-Shaped Temperature Loading 

Figures 20 and 21 show, respectively, the buckled shapes of the wing segment under actual 
temperature loading (fig. 15) with fixed and free stiffeners. Note that the buckling mode of the wing 
segment is very similar to that of the fore panel with the CL-SS condition (fig. 17-c). The calculated 
buckling temperatures of the wing segment with fixed and free supports are shown in Table 7. The 
uniform temperature loading case is also included in that table for comparison. 


Table 7. Buckling temperatures of wing segment panel 
under actual and unifonn temperature loading. 



( T ) °F 

uJcn * 

(T ) °F 

V J o/cn A 


Unifonn FX 

260 


1.00 

Actual FX 


389 

1.50 

Actual FR* 


1,324 

5.09 


* Closest to actual case 
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The wing segment model has only one row of E43 elements in the span-wise direction. Tims, the 
buckling solution termed "‘actual FR” in Table 7 may not be accurate enough to represent the buckling 
temperature of the two-dimensional wing panel. However, Table 7 can show how the “buckling 
temperature magnification factor of the first kind” 77 changes with the thermal loading and the boundary 
condition. Notice that the buckling temperature for the actual FR case is magnified by 5.09 times from 
the uniform FX case. 


The unit panel was identified as the potential thermal buckling initiation zone based on the results of 
the uniform temperature loading buckling analysis (fig. 18 and Table 6 ). Therefore, thermal buckling 
analysis of the Hyper-X wing panels may be reduced to the thermal buckling analysis of the unit panel 
without going through complex modeling of the entire wing structure. Figures 22 and 23 show 
respectively the buckled shapes of the unit panel under dome-shaped temperature loading 
(T S /T o = 0.54, Hyper-X case) with fixed and free supports. For the fixed supports (fig. 22), the unit 
panel buckled into single half-waves in two different directions, one along the radial stiffeners and the 
other transverse to these radial stiffeners. But for the free support case (fig. 23), it buckled into two 
half-waves along the radial stiffeners. Table 8 shows buckling temperatures of the unit panel under 
dome- shaped temperature loading (T S /T 0 = 0.54, Hyper-X case) with different support conditions. 
The values of the 77 are also listed in the table. Data given in Table 6 for the uniform temperature loading 
case are also listed for comparison. 


Table 8 . Buckling temperatures of unit panel; T S /T Q = 0.54 (Hyper-X case). 



Uniform FX 


Dome-shaped FX 

Dome-shaped FR 

Support 

(T 0 ) cn °F 

V 

<Jo)cn °F 

i) 

(T 0 )cn °F 

V 

SS-SS 

97 

LOO 

132 

1.36 

1.389* 

14.32 

SS-CL 

239 

LOO 

357 

1.49 

2,290* 

9.58 

CL-SS 

114 

1.00 

147 

1.29 

1,411 

12.38 

CL-CL 

295 

1.00 

412 

1.40 

2,298* 

7.79 


* Closest to Hyper-X case 

• Approaching melting range 2,375 °F ~ 2,500 °F of Haynes 230 alloy 


Note that for the fixed supports, by changing the temperature loading from uniform to dome-shaped, the 
values of 77 are in the range of 1.29 < 77 < 1.49 for all the support fixed conditions (low magnification in 
thermal buckling temperatures). However, for the dome- shaped temperature loading with free supports, 
the 77 values reached as high as 7.79 < 77 < 14.32 range, implying a great magnification of buckling 
temperatures. The dome-shaped temperature loading with the SS-SS free support condition (closest to the 
Hyper-X case) has the highest 77 value (77 = 14.32) and the CL-CL case has the lowest 77 value (77 = 7.79). 
For the SS-CL and CL-CL cases, the buckling temperatures (indicated with solid dot in Table 7) 
approached the melting range 2,375 °F ~ 2,500°F of the wing panel material. The buckling temperature 
1,389 °F of the unit panel under the dome-shaped- free-SS-SS support condition (quite close to the 
buckling temperature 1,324 °F of the wing segment panel), could be a reasonable approximation to the 
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actual buckling temperature of the aft wing panel of the Hyper-X. The buckling temperature 1,389 °F of 
the unit panel is 1.20 times the calculated peak temperature 1,153 °F of the wing panels (fig. 14). Thus, it 
is unlikely that the Hyper-X wing panels will buckle during the nominal Mach 7 flight trajectory. 

Table 9 andlO, respectively, show the buckling temperatures of unit panels under simply supported 
and clamped freely expandable boundaries subjected to dome-shaped temperature loading with different 
heat sink temperatures. Tire i] and q are also listed. 


Table 9. Buckling temperatures of unit panel under dome-shaped temperature loading; 
simply supported free boundaries. 


T/T 0 0 . 

0.1 

0.2 

0.3 

0.4 

0.5 

0.54* 

0.6 0.7 

0.8 

0.9 

1.0 

(T 0 ) cn °F 826 

872 

971 

1,067 

1,189 

1,352 

1,389 

1,458 1,732 

2,507* 

4,866* 

oo 

o 

q 

t-H 

1.06 

1.18 

1.29 

1.44 

1.64 

1.68 

1.77 2.10 

3.04 

5.89 

oo 

77 8.52 

8.99 10.01 

11.00 

12.26 

13.97 

14.32 

15.03 17.86 

25.85 

50.16 

oo 

•Exceeded incipient melting point 2,345 °F of Haynes 230 alloy. 
* Hyper-X 






Table 10. Buckling temperatures of unit panel under dome-shaped temperature loading; 
clamped free boundaries. 


T/T 0 0. 

0.1 

0.2 

0.3 

0.4 

0.5 

0.54* 

0.6 0.7 

0.8 

0.9 

1.0 

(T 0 ) cn °F 1,294 

1,385 

1,408 

1525 

1661 

1,857 

1,942 

2,141 2,763* 

4,333* 

8,447* 

oo 

o 

q 

t-H 

1.07 

1.09 

1.18 

1.29 

1.44 

1.50 

1.65 2.14 

3.35 

6.53 

oo 

77 4.39 

4.69 

4.77 

5.17 

5.63 

6.29 

6.58 

7.26 9.37 

14.69 

28.63 

oo 


•Exceeded incipient melting point 2,345 °F of Haynes 230 alloy. 
* Hyper-X 


The solid dot symbols in both Tables 8 and 10 indicate that the buckling temperatures have exceeded 
the melting range 2,375 °F ~ 2,500 °F of the wing panel material. At T S /T Q = 1 for which the 
dome-shaped temperature distribution degenerates into uniform (flat) temperature distribution, the 
buckling temperatures reached infinity, implying that thermal buckling does not occur. Note that at 
T s /T o = 0 in Tables 9 and 10 (no heat sink thermal expansion), the buckling temperatures 
{( T 0 ) cr = 826 °F, (T o ) cr = 1,294 °F } for the free simply supported and free clamped cases are much 
higher than the corresponding buckling temperatures { (T 0 ) cr = 132 °F, (T o ) cr = 412 °F} for the 
fixed support cases (Table 8 ). This results from the fact that for the free support cases, the unheated heat 
sink still can deform elastically because of the thermal expansion of the unit panel, causing the buckling 
temperatures to be much higher. 
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For the simply supported case (Table 9), T) is related to £ through 


V = 


( T ) 

v o'er 

(T ) 

v u'er 


826 « 
97 ? 


8.52^ 


( 7 ) 


And, for the clamped case (Table 10), q is related to q through 


(^) 

( r «). 


1294 g. 
295 * 


4.39^ 


( 8 ) 


The simply supported case has the numerical coefficient for q roughly twice of that for the clamped 
case (right-hand sides of equations (7) and ( 8 )). 

The data shown in Tables 9 and 10 are plotted respectively in figures 24 and 25 for easy visualization 
of how the buckling magnification factors t] and q increase with the increasing heat sink temperature 
T/T q . The TJ curve for the clamped free boundary case lies much below that for the simply supported 
free boundary case (fig. 24). Tire q curves (fig. 25) for the simply supported and clamped free boundary 
cases are quite close (Tables 9, 10). At T ^/T Q = 1 , the dome-shaped temperature distribution 
degenerates into uniform temperature distribution, and the buckling temperatures reach infinity (i.e., no 
buckling). Both figures 24 and 25 are very useful for finding the value of 77 or <q associated with heat sink 
temperature T S /T Q at certain time steps after the launching. 


CONCLUSIONS 


Heat transfer, thermal stress, and thermal buckling analyses were performed on the Hyper-X whig 
structure for the Mach 7 mission. The key results of the analyses follow. 

1. The maximum temperature differential between the upper and lower skins of the Hyper-X wing 
occurs at 89 sec after air launching from the B-52 aircraft. The profile of temperature distribution 
over the wing panel of each bay is dome-shaped. 

2. The maximum shear stress T max at the welded site is estimated as x max = 38,564 lb/in 1 2 3 4 5 , and is 
about 71 percent of the shear failure stress of the welding material. 

3. By identifying the unit panel region as the potential thermal buckling initiation zone, thermal 
buckling analysis of the Hyper-X wing panels may be reduced to the thermal buckling analysis of 
the unit panel without going through complex modeling of the entire whig structure. 

4. The estimated buckling temperature of the unit panel is (T 0 ) cr = 1,389 °F which is 1.20 times 
higher than the peak temperature 1,153 °F experienced by the wing panels under Mach 7.5 
heating. This unit panel buckling solution may be considered as a reasonable estimate for the 
buckling temperature of the Hyper-X wing panels. This result suggests that Hyper-X wing panels 
are unlikely to buckle during a Mach 7 mission. 

5. “Buckling temperature magnification factor of the first kind” r\ was established to scale up the 
panel buckling temperature in the uniform temperature loading case to estimate the panel buckling 
temperature of the dome- shaped temperature loading case. 
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6. “Buckling temperature magnification factor of the second kind” E, was established to scale up the 
panel buckling temperature in the dome-shaped temperature loading case with unheated boundary 
heat sinks to generate a solution for the same case when boundary heat sinks are heated. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, April 12, 2001 
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APPENDIX 


The temperature-dependent physical properties of Haynes 230 used for input to the thermal and 
structural models are listed in the following tables (Haynes International, Inc., Kokomo, Indiana). 


Table A-l. Thermal properties of Haynes 230. 


T op 

c p , Btu/lb °F 

k, Btu/hr-ft °F 

£ 

70 

0.095 

5.167 

0.85 

200 

0.099 

5.917 

0.85 

400 

0.104 

7.250 

0.85 

600 

0.108 

8.500 

0.85 

800 

0.112 

9.833 

0.85 

1000 

0.112 

11.083 

0.85 

1200 

0.134 

12.333 

0.85 

1400 

0.140 

13.667 

0.85 

1600 

0.145 

14.917. 

0.85 

1800 

0.147 

16.250 

0.85 
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Table A-2. 

Material properties of Haynes 230; p 

= 0.324 lb/in 3 . 


T o F 

E x 10 5 6 , lb/in 2 

<7 t x 10 3 , lb/in 2 

rt Y x 10 3 , lb/in 2 

a x 10 6 , in/in- °F 

* 

V 

70 

30.6 

125.4 

57.4 

7.0 

0.310 

200 

30.1 

(122.2) 

(55.0) 

7.1 

0.311 

400 

29.3 

(117.3) 

(51.3) 

7.2 

0.315 

600 

28.3 

(112.3) 

(47.7) 

7.4 

0.318 

800 

27.3 

(107.4) 

(44.0) 

7.6 

0.321 

1000 

26.4 

102.5 

40.3 

7.9 

0.324 

1200 

25.3 

97.7 

39.5 

8.1 

0.330 

1400 

24.1 

87.7 

42.5 

8.3 

0.332 

1600 

23.1 

63.1 

37.3 

8.6 

0.334 

1800 

21.9 

35.2 

21.1 

8.9 

0.340 

2000 

20.7* 

19.5 

10.8 

9.2* 

0.343 


* Estimated; ((Interpolated 
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Figure 1 . Hyper-X hypersonic flight research vehicle mated to the 
launch vehicle Pegasus® booster rocket. 
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Figure 2. Hyper-X hypersonic flight research vehicle flight trajectory. 
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Figure 4. Finite-difference thermal model for Hyper-X wing segment; 
welded site dimensions in inches, full contact and partial contact. 



(a) Combined case. 



(b) Separate case. 
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Figure 5. SPAR finite-element models for computations of internal radiation view 
factors. 
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(a) Altitude. 
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(b) Angle of attack. 
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(c) Mach number. 



Figure 6. Nominal Mach 7 mission flight trajectory of the Pegasus booster rocket carrying the Hyper-X 
vehicle for heating calculations. 
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Figure 7. SPAR finite-element structural model for Hyper-X wing segment. 
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(a) Uniform temperature loading. 


T 



(b) Dome-shaped temperature loading. 


Figure 9. Two types of temperature loading for the Hyper-X wing unit panel. 
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Figure 10. Nodal temperature time histories at the middle spar upper welded 
site; nodes 13, 70. 



Figure 11. Nodal temperature time histories at the middle spar lower welded 
site; nodes 53, 74. 
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Figure 12. Nodal temperature time histories in the vicinity of the middle 
spar upper welded site. 



Figure 13. Upper skin nodal temperature time histories near the middle spar 
upper welded site. 
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Figure 14. Nodal temperature time histories for node 7 of bay 1 upper skin 
and node 46 of bay 1 lower skin. 
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Figure 15. Chord- wise distributiuon of Hyper-X wing segment skin 
temperatures; t = 89 sec, 0.04 in. web contact. 
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Figure 16. Deformed shape of welded site with element shear 
stresses indicated, lb/in 2 ; Mach 7.5 heating (nominal Mach 7 
mission). 
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(Tu)cr = 167 °F 
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(a) SS-SS condition. 



(b) SS-CL condition. 



(c) CL-SS condition. 



(d) CL-CL condition. 


Figure 17. Buckled shapes of fore panel under different fixed support 
conditions, uniform temperature loading. 
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(a) SS-SS condition. 



(b) SS-CL condition. 



(c) CL-SS condition. 



(d) CL-CL condition. 

Figure 18. Buckled shapes of aft panel under different fixed support 
conditions, uniform temperature loading. 
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(Tu)cr - ^ 


(a) SS-SS condition. 



(b) SS-CL condition. 



(c) CL-SS condition. 



(d) CL-CL condition. 


Figure 19. Buckled shapes of unit panel under different fixed support 
conditions, uniform temperature loading. 
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Figure 20. Buckled shape of wing segment under actual temperature loading 
at t = 89 sec, fixed stiffeners. 



Figure 21. Buckled shape of wing segment under actual temperature loading at 
t = 89 sec, free stiffeners. 
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T 0 cr=132 °F 


(a) SS-SS condition. 


(T 0 )cr = 357 °F 




(b) SS-CL condition. 
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(T 0 )cr = 147 °F 


(c) CL-SS condition. 

010170 

_ 

(To)cr = 412 °F 


(d) CL-CL condition. 


Figure 22. Buckled shapes of unit panel under different fixed support 
conditions, dome-shaped temperature loading, TJT 0 = 0.54 . 




(a) SS-SS condition. 



(b) SS-CL condition. 



(c) CL-SS condition. 



(d) CL-CL condition. 
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Figure 23. Buckled shapes of unit panel under different free support 
conditions, dome-shaped temperature loading, TJT a = 0.54 . 



Figure 24. Plots of buckling temperature magnification 
factors of the first kind, t ] , as functions of boundary heat 
sink temperature TJT 0 . 



T s /T 0 010177 


Figure 25. Plots of buckling temperature 
magnification factors of the second kind , g , as 
functions of boundary heat sink temperature TJT a . 
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